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Abstract 

We study the previously constructed Riemann problem whose so- 
lutions correspond to equilibrium configurations of black holes. We 
evaluate the metric coefficients at the symmetry axis and the interac- 
tion force between the black holes. 



1 Introduction 



In the previous work we studied the axially symmetric stationary vac- 
uum solution of Einstein equations describing the equilibrium configuration 
of rotating black holes. The equilibrium configuration was understood to 
mean a stationary solution possessing a disconnected event horizon. All such 
solutions satisfy a certain boundary-value problem for a system of elliptic 
nonlinear equations that can be easily obtained from H , where the regularity 
conditions for the symmetry axis and the event horizon were first formulated. 
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Our investigation showed that this boundary-value problem reduces to a ma- 
trix Riemann problem with a rational conjugation matrix. We also noted 
that solutions corresponding to stationary equilibrium configuration of black 
holes are likely to develop a conical singularity on the symmetry axis, which 
is the reason why the black hole do not fall on each other. We note that for 
rotating black holes, the presence of conical singularity for arbitrary values 
of the solution parameters is not quite obvious from the physical standpoint, 
because it can be assumed that the rotation of black holes leads to repulsion 
forces that can compensate the gravitational attraction. 

The conical singularity itself allows us to introduce the notion of the 
interaction force between the black holes. We consider this in more detail in 
the next section. The analytic expression for the interaction force between 
Schwarzschild black holes has long been known ||. The main result of the 
present work is our generalization of this to the case of rotating black holes. 



2 The boundary-value problem 

In cylindrical coordinates, the axially symmetric metric is 

ds 2 = -Vdt 2 + 2Wdtdcj) + Xd<p 2 + ^e\dp 2 + dz 2 ), 

P 2 

where the metric coefficient depend only on p and z. The Einstein equations 
are then divided into two system of nonlinear equations, the first of which 
can be brought to the form 

{pg, P g~ l ), P + {pg,zg~ x ), z = o, g = ( ^ x ) , det# = -p 2 , (2.1) 



and the second one to the form 



(Hj 2 e P )),C = Y p ( 1+ ^(PHg- 1 ?) , d ( = \{d z - id p ). (2.2) 

Let the event horizon have N disconnected components, let Zi, . . . z 2 n 
be the ^-coordinates of the intersection points of the event horizon and the 
symmetry axis, and let fli be the angular velocity of the ith black hole. In 
the neighborhood of the ith black hole, we choose the coordinate system 

2 / \2 2wi 2\ z 2i ~ z 2i-l 

p = (A -mj(l-/i), rrii = , 
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z = A/i, \n\ < 1, 

in which the regularity condition of the symmetry axis and the event horizon 
can be written as 01 



1 n 4 \ /l ON _ /(A - mf)U(A,p) pW(A,p) 



iy y Va 1/ V p W(X, p) (l-p 2 )X(A,p) 



(2-3) 



where X and are smooth functions not equal to zero. Using (2.3), we can 
easily obtain the boundary conditions for the system (2.1) 

P9, P g~ l = ([J ^), p^o, zer, (2.4a) 

CUpg.pg-'Clr 1 = [ 2 {1) °),p-0«G/ <i n i =(J (2.46) 

p^- 1 = 0(l), p^O, zGi?. (2.4c) 
where T is the symmetry axis consisting of N + 1 connected components, 

r = |Jr i = i2\U^ U = i,...,N + i,i = i,...,N), h = (z 2l _ u z 2i ). 

The symbol 0(1) denotes uniformly bounded functions on the corresponding 
interval. We impose the condition at infinity 
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p 2 0(- 3 ), X = p 2 (l + 0(-)), r = + Z 2. (2.5) 



Boundary conditions (2.4) and (2.5) are sufficient to construct solutions of 
the first group of Einstein equations (2.1). 

We now consider the properties of solutions of system (2.2). Using (2.3), 
we can easily verify that for p = and z G T, we have 

d z (3 = 0, /3\ r . = bi, (2.6) 

where hi are some constants. The condition for the solution to be regular on 
the symmetry axis { X \x' a ~ > ^ on ^ ne symmetry axis) is satisfied if and only 
if Eq.(2.3) is supplemented by the additional requirement that 

h = 0. (2.7) 
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For asymptotically flat solutions (those for which (3 — > as r — > oo) condition 
(2.7) is automatically satisfied on the first and the last components of the 
symmetry axis but cannot be satisfied on the remaining components (see 
[0- H H' 0)- I n ^ ne general case, the last statement is not yet proved. In this 
work, we evaluate the constant hi, but the corresponding analytic expression 
is unfortunately too complicated to analyze the solvability of Eq.(2.7). 
The hi parameters also have their own interpretation: the quantity 



can be considered as the interaction force between the black holes. Here and 
what follows, we assume that h\ = &jv+i = 0. We now comment on the origin 
of this interpretation. If condition (2.7) is not satisfied, the curvature tensor 
becomes singular in the neighborhood of the corresponding symmetry axis 
components. We assume that the singular points are filled with "matter", 
which is precisely what prevents the black hole falling on each other. The 
tension in z direction per unit surface is T(e z , e z ) where T is the energy- 
momentum tensor and e z is the normalized vector parallel to d/dz. It is 
natural to define the interaction force of the black holes as the integral 



where S £ is the 2-surface coordinatized by p0(with p < e) and ds is the 
corresponding area element. Evaluating this integral, we obtain Eq. (2.8) 



To conclude this section, we give a brief derivation of representation for 
hi We bring Eq.(2.2) to the form 



(2.8) 




I 




X,2 , {d p Yf-{d z Yf 



X 2 



) 



(2.9a) 



d z (3 = d z InX (pd p In X - 2) + 



pd z Yd p Y 
X 2 



(2.96) 



where 



dY 



1 



*{XdW -WdX), (*dp = dz,*dz 



dp). 



P 
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Let C £ be a curve connecting T i+ i and Tf 

p — esmr, z = — y e 2 + mf cost, < t < 7r. 

Then 

k +1 -h= [ d{3. (2.10) 

The left-hand side of (2.10) is independent of e. Taking e to zero and using 
(2.3) and (2.9), we obtain that 

bi+i -bi = -2 (\nX(z 2l ) - \nX(z 2l ^)) , (2.11) 

where X(z 2i ) = X(rrii, 1) and X(z 2 j-i) = X(jrii, —1). In the next section, we 
use Eq.(2.11) to evaluate 6j. 

3 The Riemann problem 

We showed in [Q that for every solution of boundary- value problem (2.4), 
(2.5) there exists a piecewise analytic matrix x( a; ) satisfying the conjugation 
condition on the imaginary axis 

X _(uO = X + (")(J l)T(k)(l yj, k = z + (^-p')/2uj (3.1a) 

and normalized at infinity by 

X(w) -+I, uj -> oo. (3.16) 
The rational matrix T(k) is defined as 

T(k)=D N+1 T N D N ...T 1 D 1 , (3.2) 

where 

t _ 2M^_ _ 4M o.\ /x n . \ 

Z}= (J f'V (3.3) 

V (k-z 2 j)(k-Z2j-\) k-z 2 j ' 

The Mj Lj parameters have the physical interpretation of the respective full 
mass and angular momentum of jth black hole and are related by 

Mj = rrij + 2QjLj, 2m j = z 2 j — z 2 j-i- (3-4) 
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An additional requirement that the matrix T{k) be symmetric leads to 2N+1 
nonlinear algebraic equations on the parameters Dj, Lj and Qj, which deter- 
mine Dj and Lj as functions of Qj and Zj. In particular, 

N+l N 

X n, • y.m.Uj o. (3.5) 

3=1 3=1 

Let d\ — Di, d i+ i = di + D i+ i + 4Mjf2j and 

t 2 ^(k) = C-j, ( I t 2j (k) = C^( 1 ^^)C^ 

\2njik-z2j-i) L J \u i / 

where 

1 -cL\ ( o a- 



Then applying (3.5) and easily verified identity 

1 -dj \ -d,- - 4M^- 



T Ho i J - ^ n 

we derive that 

T(k) = t 2N {k)t 2N ^{k) . . . hik^ik). (3.6) 
We now note that the matrix 

has no singularity at uj = and tends to identity matrix as iq — > oo, 



1 0\ . , IX /l 



o w /iw(*)*y-i(*) Vo 1/w 



1 -8Mj(2MjQj + dj)/p 



^=o V 1 

(3.7) 

Therefore, the same properties are shared by the matrix 

and Riemann problem (3.1) is correctly posed for any Dj,Lj,Vtj satisfying 
(3-5). 
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The only singularities of the conjugation matrix are simple poles at points 
,± 



= ( Zi -z)± y/(zi- z) 2 + p 2 , Ui = -P 2 /loI 



It now follows from Eq.(3.1) that X±{ UJ ) are rational functions with simple 
poles at the points uf (uf > 0, uj~ < 0), 

2N A± 

x± (a;) = J + £ (3.8) 

3=1 U - 

with Af being independent of uo. It follows from the unimodularity of T(k) 
that detx±(o;) = 1, and det Af = 0. Let £j be the eigenvector of Aj; then 

A7fc = 0, Aj = aj %, C; ( ij.ij). (3.9) 

with aj being a column vector. We further note that is unimodular if 

and only if there exists a constant V{ such that 

en = Vi Iim_x-M^i, (3.10) 
Equation (3.9) and (3.10) yield a system of linear equations for a^, 

a% = Viii + Vi V _ a fc . (3. 10a) 

The parameters t> j and £j can be easily determined from conjugation matrix. 
We define the matrix Gj(fc) = U-i(k) . . .t\(k) and the vectors 

= C 2j (q ) ' C2 i-! = ^-l ( J ) • 

6= (J j^Gr^*)*, (3.11) 



Then 



and also 



(3.12) 
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where 



,o __ „,o 1 1 



v 2j ~ + > V 2j-1 ( > - 

^27 — ^27 & y W 27-l — U 2 



J 2j~^2j «j - <"2i-l 

and 

Ui = 1 + ( w r - ^^G^Gr 1 ^. 
We also note that 

£f = u;r<G^)( ly ,°- 

The approach adopted here to solve the Riemann problem is similar to 
one used in [Q . As is clear from what follows, to construct an exact solution 
to nonlinear system (2.1) we can take T to be any symmetric matrix. This 
method for deriving exact solutions is evidently different from the method 
used in P, fTofl; however, it does not lead to new solutions of (2.1). In par- 
ticular, the solution that we investigate here is in the class of 2A r -solitons in 
Minkowski space background |l(J . 

We now show how the solution of (2.1) can be reconstructed from the 
solution of Riemann problem (3.1). Let 

Di = d z 9 2U; 9 <9^, D 2 = d p + <1Up - d u . 

uj z + p z uj z + p z 

Because D\k = D 2 k = 0, it follows from (3.1) that "logarithmic" derivatives 
of the functions 



*±(") = x±H(J °) 



have no other singularities in addition to the simple poles at ±ip. Therefore, 

D 1 *. M r'M - zi£M^mM + ¥^ti£m-i£> (3 . 13a ) 

i>,<M-0*-(«0 = "Wri^M + AM-^-h'Hp) (3m) 

v ' v ' uo-ip u + ip v ' 

We now take into account that T(u) = T{u) and define the real matrix 

g = -x-(o)(l _° p2 )- (3-14) 
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From (3.13) with oj = 0, one obtains 

d c gg~ l = dJi!.{ip)^Z\ip), B c gg- 1 = d^.{-ip)^!Z\-ip). (3.15) 

The compatibility condition for Eqs.(3.13) is given by system (2.1), and there- 
fore, Eq.(3.14) is a solution of (2.1). 

From here on, we assume that the parameters dj and Mj are determined 
from the requirement that T(k) be symmetric. Then 

N 

Mi{di + = 0, ( lim k(T 21 (k) - T 12 (k)) = 0) (3.16) 

i=l 

and X-(0) = X+(0) ( see (3-7)). Using the last identity and the uniqueness of 
the solution to the Riemann problem, we obtain the reduced problem 

^H = - 9 x; 1 (7M(J _ 1 ° /p2 ) (t(-p» = th), (3.17) 

where T denotes the transposed matrix. Taking u either to zero or to infinity 
in (3.17), we see that g is a symmetric matrix. 

We now investigate the properties of solutions of the Riemann problem 
as p — > 0. Let z G T m+ i, then 

c^-> 2{zi-z), i>2m + l, p^O (3.18a) 

->2(zi- z), i<2m, p^O (3.186) 

Using the formula uof = —p 2 /uoZ, we can easily verify that all the coefficients 
of linear system (3.10) also well defined for p — 0, and it is therefore natural 
to assume the existence of the limit lim p ^o a i- This assumption ensures the 
existence of all the limits used in what follows. 
Introduce the matrix 




Q l (k)=t 2N (k)...t t (k), T(k) = Qi(k)Gi(k). (3.20) 
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We can see from (3.18) that 



limxiM = /+ -. (3.21) 



Further, let 

In view of the identities 



1 0\~ ~ (I 



1 0\~ - (I 



1 %M j (d j + 2M j tt j ) 
1 



we have 



X + (0) = X 2 (0), X 2 (co)=(J «S*>"iW + ™A>). (3. 22) 
In terms of Xi an d X2, Riemann problem (3.1) can be rewritten as 

X 2 H = XiH(J °)Q 2 ^ + iG 2m+1 (J (3.23) 



and reduced problem (3.17) as 



VU VP / i=l ^ — i=2m- 



1 \ ^ Bf ™ B 



J i i=2m+l W-Ui 



(3.24) 

Taking p — > in (3.23) and taking into account (3.21), we see that 



(3.25) 



where 



7m+iM = Q^ +1 (z + co/2)G 2m+1 (,2 + w/2). (3.26) 
The function X2(^) is nonsingular at u — 0. This is possible only if 



Vo 
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Finally, 

v (Ci\ = ( Vt^iW (7^ 2 + i(0R 7 ^ +1 (0) - 7^ 2 + i(0)^ 7 ^ +1 (0)) h% + M \ 
XA) { 7 - +1 (0) J 

(3.26) 

Recall that x_(0) — X+(0) — X2(0). We now note that 7 m+ i(0) is a symmetric 
real matrix with the unit determinant and also that 7„ +1 (0) — > 1 as z — > oo; 
therefore, 

7^ 2 +1 (0)>0. (3.27) 

Since 



A 2 ^ lr - = *T3^lr< = ^(M, 



and A 2 — mf — > (z — 2 2 i)(z — z 2 j_i) as p — > 0, we find 

X(z 2i ) = mires^-Tf+^O, 2), X(^-i) = -m i res aa ._ 1 'yf 2 (0, z). (3.28) 

It follows from (3.27) that 

res^7, 2 i(0, 2) > 0, res^_ l7 f (0, z) < 0. 

Equations (2.11) and (3.28) give the sought-for analytic expression for the 
constants hi and the interaction force (2.8) of black holes. 

4 Interaction force between two black holes 

For two black holes, the symmetry axis has three connected components. We 
normalize f3 such that b± = 63 = 0, then 

F= U*M-l), (4.1) 

and 

X(z 2 ) = m 1 res Z2 7 2 22 (0, z), X( Zl ) = -m 1 res Zl7l 22 (0, (4.2) 

Expression (4.1) depends on di,fli, L iy Mi and Zj. These parameters are not 
independent and must satisfy the system of nonlinear algebraic equations 

res Zi Ti 2 (A;) = res Zi T 2i (/c), M; = m, + 2fi;Lj. (4.3) 
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It can be easily seen from (3.16) that 

di = -2Mifii + M 2 d, d 2 = -2M 2 Vl 2 - M x d. (4.4) 

Using (4.3) and (4.4), we can eliminate parameters di,d 2 ,Qi and fl 2 from 
the expression for the interaction force (4.1). Unfortunately, the formula 
thus obtained is extremely cumbersome. The situation is greatly simplified 
in one particular case, however. Namely, let the black holes rotate in opposite 
directions and have the same total and irreducible masses: 



We set 



Then 



Mi = M 2 = M, mi — m 2 = ttl, Qi = —Q 2 = Q. 

Zi = —mi, z 2 = m i, z 3 = R — m 2 , z 4 = R + m 2 . (4.5) 



R -+- 2M 

L = 2M 2 (m + M) VI, M = m + 2m, L x = -L 2 = L, d = (4.6) 

R — 2M 

and 

M 2 

F = . (4.7) 

By definition, R > 2m. Using (4.6), we can easily show that 2M < R (where 
the equality takes place only in the limiting case of R = 2m); therefore , F > 
and tends to infinity as R — > 2m. We also observe that the interaction force 
in the case under consideration is exactly the same as for two Schwarzschild 
black holes with identical masses (equal to M). The interaction force of two 
nonrotating black holes (Hi = £l 2 = 0) is given by 

F = ^ . (4.8) 

R 2 - {mi + m 2 ) 2 1 J 

It is easy to see that in the cases we have considered, F tends to the 
Newtonian limit as R — > oo. This result can be generalized; the estimate 

'-^ + °£> <-> 

holds as R — > oo. To derive (4.9), it suffices to use the approximate solution 
of (4.3): 

2L 2 Mi + 2LiM 2 1 / 2L 1 Mf + 2L 2 M 1 2 

d — — -|- 4 i7i -|- Lio -\- 

MiM 2 R 2 V M X M 2 
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a L '" ! 



2Mf(M i + m i ) 
and 

1 4M 2 8M| 1 . 1 4Mi 8M X 2 ^. 1 . 

We now discuss the general properties of constraints (4.3) on the angular 
velocities and angular momenta. We especially note that the equality L 2 = 
Q2 = cannot be satisfied for any values of the remaining parameters with 
L\ 7^ 0. This means that if L 2 = and the second black hole does not 
contribute to the total angular momentum, an approaching observer would 
still see this black hole rotating with a certain angular velocity. We consider 
this behavior quite natural. Indeed, the notion of the angular velocity is 
related to the behavior of geodesies in the neighborhood of black hole: every 
test particle approaching the black hole is involved in its rotation, and the 
angular velocity of this rotation becomes equal to the angular velocity of the 
black hole when the particle reaches the event horizon. That the angular 
velocity is nonzero is therefore determined by the fact that the first black 
hole involves the second one in its rotation. This explanation corresponds to 
the behavior of f2 2 at large distances: if L 2 = then Q 2 — 0(1/ R 3 ). 

We return to the exact solution (4.6). The foregoing discussion allows us 
to take the angular momenta of each black hole as the independent param- 
eters. We then have Q — > and M — > m as R — > 2m which means that as 
black holes approach each other, they loose the full mass and "slow down" 
each others rotation. 

Evidently, the realistic solution for the configuration of two black holes 
cannot be stationary. As we have seen however, the stationary solutions 
admits a sufficiently realistic physical interpretation and demonstrate, at 
least qualitatively, the effects caused by the interaction of black holes. In 
our opinion, this is because the Einstein equations determine not only the 
gravitational field but also the law of motion of material bodies [JTT] . 

This work was supported by RFBR grant No. 96-01-00548. 
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